Decomposition of entanglement entropy in lattice gauge theory 
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We consider entanglement entropy between regions of space in lattice gauge theory. The Hilbert 
space corresponding to a region of space includes edge states that transform nontrivially under 
gauge transformations. By decomposing the edge states in irreducible representations of the gauge 
group, the entropy of an arbitrary state is expressed as the sum of three positive terms: a term 
associated with the classical Shannon entropy of the distribution of boundary representations, a term 
that appears only for non-Abelian gauge theories and depends on the dimension of the boundary 
representations, and a term representing nonlocal correlations. The first two terms are the entropy 
of the edge states, and depend only on observables measurable at the boundary. These results are 
applied to several examples of lattice gauge theory states, including the ground state in the strong 
coupling expansion of Kogut and Susskind. In all these examples we find that the entropy of the 
edge states is the dominant contribution to the entanglement entropy. 
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I. INTRODUCTION 

In quantum field theory, regions of space are subsys- 
tems, and the entropy of these subsystems, virhich for pure 
states is the entanglement entropy, gives important infor- 
mation about the quantum state. The entanglement en- 
tropy plays an important role in black hole physics, where 
it governs the one-loop corrections to the Bekenstein- 
Hawking entropy, and mayin fact be responsible for the 
entire black hole entropy 

In addition to its role in black hole statistical mechan- 
ics, entanglement entropy has been found to play a role 
in the study of phase transitions. Entanglement entropy 
can be used to test for the presence of long-range order, 
even when an order parameter is not known [5|, |6| and 
can detect the presence of topological phases 0, Q ■ In 
gauge theories, entanglement entropy may be useful in 
studying the deconfining phase transition [l^ . 

In gauge theories, the entanglement entropy is compli- 
cated by the fact that states are not precisely localizable 
in space. The result is the the Hilbert space correspond- 
ing to a region of space includes edge states that con- 
tribute to the entanglement entropy US] • These edge 
states are similar to the "would-be pure gauge" de gree s 
of freedom in (2-|-l)-dimensional quantum gravity [iSj. 
Recall that in 2-1-1 gravity there are no local degrees of 
freedom in the bulk, yet the usual thermodynamic argu- 
ments suggest that black holes have an entropy propor- 
tional to the length of the horizon. In the approach of 
Ref. |1^, the horizon has local degrees of freedom and 
it is these degrees of freedom whose entropy is given by 
the Bekenstein-Hawking area law. While the precise re- 
lation between these two notions of boundary states is 
not clear, their similarity suggests a relation between the 
entanglement entropy and the boundary state counting 
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method of Ref. 

It is well-known that the entanglement entropy is ul- 
traviolet divergent, so a regulator is needed in its defi- 
nition. In the continuum, standard methods for calcu- 
lating the entanglement entropy at one-loop order fail 
for gauge fields ^^y. The Hamiltonian method does not 
work because the eigenfunctions of the boost generator in 
Rindler space are not square-integrable. The Euclidean 
conical deficit angle method leads to a "contact term" 
that makes the result negative and therefore not identi- 
fiable with entanglement entropy, which is a manifestly 
positive quantity. In this paper we consider the entangle- 
ment entropy for lattice gauge theories. This allows us to 
regulate the entanglement entropy while clearly exposing 
the role of gauge invariance. 

Entanglement entropy has been considered before for 
certain classes of states in lattice gauge theory [l^, ITB - 
[Tsj . This paper adopts the Hamiltonian formulation of 
lattice gauge theory, rather than the replica method fl^ 
that has typically been used in numerical calculations 
of the entanglement entropy. The replica method re- 
lates the entanglement entropy of the ground state of 
a given theory to the partition function computed on an 
n-sheeted cover of the Euclidean spacetime. Our results 
do not use the replica method, so they do not require 
the state to be expressed as a Euclidean path integral, 
though we agree with results obtained using the replica 
method where the latter is applicable. 

Closely related to lattice gauge theory is loop quan- 
tum gravity, which is formulated as an SU(2) lattice 
gauge theory on a superposition of lattices. Although 
this paper will not discuss loop quantum gravity, entan- 
glement entropy in loop quantum gravity was discussed in 
Refs. l21j, and we expect the techniques of this paper 
to generalize easily to a superposition of lattices. We note 
also that the Hilbert space of edge states in SU(2) lattice 
gauge theory is closely related to the Hilbert space of the 
SU(2) Chern-Simons theory whose states are counted in 
the loop quantum gravity derivation of black hole entropy 
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[2^1231. 

We now briefly summarize our result. Consider a lat- 
tice whose set of nodes N is divided into two disjoint 
sets A and B whose union is all of N. In a lattice the- 
ory where the degrees of freedom live on the nodes, the 
Hilbert space associated with a set of nodes is simply the 
tensor product of the Hilbert spaces of each individual 
node. This leads to a tensor product decomposition of 
the whole Hilbert space as H = Ha (^Hb- In a lattice 
gauge theory, the degrees of freedom live on the links, so 
there is not such a simple tensor product decomposition. 
However, following Ref. we can define a Hilbert space 
Ha by splitting the links that cross the boundary. Along 
each link I with one end point in A and one end point in 
B, we insert a new vertex on the boundary and divide 
the link into two smaller links, one associated with region 
A and one associated with region B. The Hilbert space 
Ha then consists of functionals of the connection on the 
links in A that are invariant under gauge transformations 
that act on the nodes in the interior of A, but not on 
the boundary nodes. In restricting the gauge symmetry, 
degrees of freedom that were previously pure gauge are 
promoted to physical degrees of freedom. The new de- 
grees of freedom are edge states that are associated with 
the boundary vertices and transform nontrivially under 
gauge transformations acting on the boundary. They are 
the lattice analogue of the continuum edge states studied 
in Ref. 

The Hilbert space H is not equal to Ha^Hb, since the 
former is invariant under all gauge transformations, and 
the latter is invariant under only those gauge transfor- 
mations that act trivially on the boundary. Thus instead 
of an isomorphism of Hilbert spaces, we have the embed- 
ding 



H Ha®Hb- 



(1) 



The entanglement entropy of any state in H can be de- 
fined by embedding the state into Ha®Hb- Letting p de- 
note a state, represented as a density matrix in Ha®Hb, 
the reduced density matrix of system A is the partial 
trace pA = tr^ (p) , and the entanglement entropy is 

"Scntanglcmont — S(^Pa) (2) 

where the function S{p) is the von Neumann entropy. 



S{p) = -tr(plnp). 



(3) 



The states in H are invariant under all gauge transfor- 
mations, including those acting on the boundary. The 
reduced density matrix pA associated with a gauge- 
invariant state is then also invariant under the group of 
boundary gauge transformations, which acts nontrivially 
on Ha- When decomposed into irreducible representa- 
tions of the group of boundary gauge transformations, 
the matrix pA takes the form of a direct sum of tensor 
products. Using properties of the von Neumann entropy 
under direct sum and tensor product, we decompose the 
entanglement entropy of a generic state in lattice gauge 
theory as a sum of three positive terms 



• the Shannon entropy of the distribution of bound- 
ary link representations found in Ref. 12], 

• the weighted average of the logarithm of the di- 
mension of the boundary representations found in 
Ref. [2?], and 

• a third term that captures nonlocal correlations of 
the bulk field. 

The first two terms are purely local to the boundary, 
and together they capture the entropy of the edge states. 
Our result directly generalizes the results of Refs. [H, 
[2ll |. We will derive the result in Sec. ITTl and give several 
applications to specific states of lattice gauge theory in 

Sec, mil 



II. 



ENTANGLEMENT ENTROPY IN LATTICE 
GAUGE THEORY 



We first review the kinematics of Hamiltonian lattice 
gauge theory and spin network states [1^. Consider 
a lattice consisting of a set of nodes and a set L of 
oriented links, and let G be a gauge group that is either 
a compact Lie group, or a discrete group. A field con- 
figuration is an assignment of group elements ui to links, 
and a gauge transformation is an assignment of group 
elements gn to nodes, which acts on ui as 



s{l) 



(4) 



where s{l) and t{l) are respectively the nodes at the 
source and target of the link I. The Hilbert space H con- 
sists of square-integrablc functionals of the holonomies 
ui that are invariant under gauge transformations. An 
orthonormal basis for H is given by a generalization of 
the spin network states [25l - l27l | . A spin network consists 
of an assignment of irreducible representations R = {ri : 
I G L} to each link, and intertwiners / = {z„ : n G N} 
to each node. Each intertwiner j„ is a G-invariant linear 
map between representation spaces 



ri 



ri 



(5) 



The spin network state \S) associated with a spin network 
S is the functional obtained by taking the representation 
ri of the group element on each link I , multiplying by 
•\/dim(r), and contracting the free indices with the inter- 
twiners in (25i] . 



{S\U)^\Q^y/dh^riiui)\ol(^in\. (6) 

\leL ) \neN ) 

The intertwiners are chosen to be orthonormal in the 
inner product, 



(^1,12) = tr(ii4) 



(7) 
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so that the resulting spin network states form an or- 
thonormal basis of H [28j] . 

We now describe the tensor product decomposition of 
the Hilbert space, which was described in Ref. l^j for 
Abehan lattice gauge theories. Let A be a region of space, 
which on the lattice will mean a subset of the nodes. The 
configuration space of Ti, consists of holonomies on all 
links of the lattice. The links can be divided into three 
sets: La is the set of links with both end points in 
Lb is the set of links with both end points in B and Lq 
is the set of links that cross the boundary. In order to 
partition the degrees of freedom of the boundary links 
between TIa and 'Hb, we split each boundary link into 
two at the boundary, such for each link in Lg there is a 
new link in Lqa and one in Lqb- The Hilbert space Ha 
is then defined as the square-integrable functions of the 
holonomies {ui : I G La U Lqa} invariant under gauge 
transformations acting at nodes in the interior of A (but 
not under gauge transformations acting on the bound- 
ary). Then for each link I g Lg, the holonomy ui can be 
obtained as a product of a holonomy in Lqa and one in 
Lqb, and we define the product map 

TT : UsA X UaB Uq. (8) 

The puUback map tt* then gives an embedding 

TT* -.U'^'Ha^Hb. (9) 

For example, if there is just one boundary link I split into 
links li and I2 such that t{li) — 5(^2): then the map tt is 
given by 

7r(ui,U2) = U2-U1, (10) 
and the pullback 7r*-0 of a function : G ^ C is given 

by 

(7r*'0)(wi,M2) = ■0(U1W2). (11) 

This embedding preserves the norm of the state, a fact 
which follows from the G-invariance and normalization 
of the Haar measure. 

To specify a spin network S, we specify all its repre- 
sentations and intcrtwincrs 

S^{RA.RB,RaJAjB). (12) 

Just as the space T-L is spanned by spin network states, 
the space T-La is spanned by open spin network states 
[29j . An open spin network Sa is specified by 

SA^{RA,RdjA,M) (13) 

where M — {mi : I S Lq} is a set of vectors in the bound- 
ary representation spaces, mi G r; if the link / points 
inward at the boundary, or in the dual representation 
mi S fi if the link I points outward (for unitary represen- 
tations, the dual representation r and complex conjugate 
representation r* coincide) . The open spin network state 



\Sa) is defined just as in Eq. ([6]), except that the ex- 
tra free indices associated with the boundary vertices are 
contracted with the vectors mi. The open spin network 
states form an orthonormal basis of T-La provided the vec- 
tors mi and intertwiners i„ are chosen to be orthonormal. 

As shown in Ref. |;21], under the embedding tt*, the 
spin network state |S') maps to 

7r*|5)= n , Y.\^a)®\Sb) (14) 

where Sa is given by (fT3|) . and Sb = {Rb, Ra, Ib, M*), 
where the vectors M* are dual (complex conjugate) to 
the vectors M . In Eq. (jl4p . mi ranges over an orthonor- 
mal basis of r; . This decomposition follows from inserting 
a resolution of unity at each point where a link crosses the 
boundary, and the factors of l/-\/dim(r;) arise to cancel 
the extra factors of dim(r;) in Eq. ^ that come from 
splitting the boundary links. 

We now consider an arbitrary gauge-invariant state \4>) 
expressed in the spin network basis, 

s 

Using the decomposition (jl4p . the reduced density matrix 
for region A is 

Ra,RBjB,M 

where S and Sa are given by p2|) and ((T3| . and S' and 
S'j^ are given by 

"5" = iR'AiRB,Rd,l'AilB), S'a = {R'A,Rd,l'AiM). 

(17) 

The sums over intertwiners in Eq. (|16|) are taken over an 
orthonormal basis of the space of intertwiners compatible 
with the representations incident on each node. In the 
case where there is no such intertwiner, the sum is zero. 

There are two features worthy of note about Eq. (|16p . 
First, the set of representations Rg is always the same 
for Sa and S*^. This means that the matrix pA has no 
off-diagonal terms that mix different boundary represen- 
tations. Second, the coefficients in Eq. are indepen- 
dent of M, so within each representation the M degrees 
of freedom are in a maximally mixed state. 

This structure of the reduced density matrix can be 
seen from group theory. The boundary gauge trans- 
formations form a group G" where n is the number of 
boundary links, and is represented unitarily on Ha- Such 
a representation R{g) can always be written as a direct 
sum of irreducible representations: 

^(5) -0K5)®lnW (18) 
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where g S G", r runs over all irreducible representations 
of G", and n{r) is the multiplicity with which the ir- 
reducible representation r appears in the representation 
R. The reduced density matrix pA comes from a gauge- 
invariant state, so it must commute with the representa- 
tion R. To commute with R{g) for all g, pA must take 
the form 



PA 



e 



n 

-^dim(r) 

dim(r) 



«) PA{r) 



(19) 



where pa^t) is a density matrix of dimension n{r). 

To see more explicitly how the density matrix de- 
composes into representations, it is useful to divide the 
Hilbert space T-La into an edge Hilbert space and a bulk 
Hilbcrt space such that the states of the boundary Hilbert 
space are labeled by M), and the bulk Hilbert space 
is labeled by {RatIa)- This decomposition of the Hilbert 
space is 



-Ra 



((g) n) 


®nA{Rd) 


. \leLa / 





(20) 



where %A{Rd) is spanned by states \Ra,Ia) that are 
compatible with the assignment of representations i?g to 
the boundary. 

In the decomposition (j20[) the open spin network states 
\Sa) and jS*^) can be written as 



\Sa) = \Rd)®\M)®\RA,lA) ■ 
\S'a) ^\Ra)®\M)®\R'^,rj,) 



(21) 
(22) 



so that their outer product takes the form 
\Sa){S'a\ = \Rd){Rd\ ® \M){M\ ® \Ra,Ia){R'aJ'a\ ■ (23) 

Substituting Eq. (|23|) into the reduced density matrix 
([in]) and rearranging terms yields 



PA 



= Y,p{R9)\R9){R9\®[Y, 



\M){M\ 



Ra 



M 



n 



dim(r;) 



^PA{Rd) 
(24) 



where p{Rd) is the probability of distribution of repre- 
sentations on the boundary. 



Ra.^RbJaJb 

and pA{Rd) is the reduced density matrix 



(25) 



PA{Rd) 



E 



IaJ'aJb 



p{Rd) 



\RaJa){R'a.I'a\- (26) 



The factor p{Rd) is included in the definition of pA{Rd) 
to maintain the unit trace condition. 

Since the same Rg appears in both the ket and the 
bra in the first tensor factor in Eq. (|24l) . the state does 



indeed lie in the direct sum Hilbert space ((20|) . Moreover, 
the second tensor factor in Eq. ([^ is proportional to the 
identity matrix, so the density matrix pA can equivalently 
be written 



PA 



^P{Ro) 



tr 



\ieLg 



dim(n) 



(EipAiRd) 



(27) 



The structure of the reduced density matrix (|27p allows 
us to simplify the entanglement entropy by using prop- 
erties of the von Neumann entropy under direct sum and 
direct product. Let p„ be positive real numbers sum- 
ming to one, and pn density matrices on Hilbert space 
Hn- The von Neumann entropy of a weighted direct sum 
is 

S (0PnP„j = Hipn) + {S{pn)) (28) 

where (•) denotes expectation value with respect to the 
probability distribution p„, and H{pn) is the Shannon 
entropy of this distribution (the classical analogue of the 
von Neumann entropy). 



(29) 



Under a tensor product, the von Neumann entropy is 
additive. 



S{pi®P2) = S{pi)+S{p2). 



(30) 



Finally, the maximally mixed state of dimension n has 
entropy Inn, 



S'(l„/n) = Inn. 



(31) 



Applying the properties of the von Neumann entropy 
JH) , (PH)) and (I5T|) to the reduced density matrix pA (^7)) 
gives the entropy as the sum of three positive terms. 



E 

l£LgA 



S{pa) = H{p{Re)) + V (Indim(n)) + {SipA{Ro))) 



(32) 

where (•) denotes expectation value with respect to the 
probability distribution p{Rg). Eq. (j32]) is the main re- 
sult of this paper. Individual terms in this expression 
have appeared before: the first term appeared in Ref. [11] 
where it was derived for a specific class of states (see sec- 
tion |IlT^, and the second term appeared in Ref. [2l| as 
the entropy of a single spin network state. 

The first two terms of Eq. ([5^ depend only on the 
distribution of the boundary representations, and in this 
sense are purely local. The second term is a sum over 
boundary links, and so is extensive on the boundary. The 
first term is not extensive, but will be approximately 
extensive as long as the correlations between different 
representations are local. The effect of correlations is 
always to decrease the entropy, so we can obtain an ex- 
tensive upper bound by neglecting these correlations. If 
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we assume that the statistics of each edge are the same 
(which would be the case for states with a discrete trans- 
lation and rotation symmetry, such as the ground state of 
a translation- and rotation-invariant Hamiltonian) then 
the upper bound depends only on the probability distri- 
bution of representation on each edge p{r) and is given 
by 

'S'boundary < n(^H{p) + (Indim(r)) ^ , (33) 

with n the number of boundary links. 

In principle either of the local terms can be larger. 
For example, in a state sharply peaked on spin net- 
works with high-dimension representations, the second 
term will dominate. In a state that is a superposition of 
many spin networks with low-dimension representations, 
the first term will dominate. In particular, an Abelian 
theory has only one-dimensional representations so the 
second term in Eq. (|32|) vanishes. 

The third term in Eq. (|32|) is the most difficult to char- 
acterize. It includes the effects of correlations between 
distantly separated degrees of freedom, and in general it 
is not bounded by the area of the boundary. However 
we will see that, for the classes of states considered in 
Sec. mil this term is either vanishing or much smaller 
than the local terms. 



III. EXAMPLES 

We now consider several examples of states whose en- 
tanglement entropy can be calculated using our method. 



A. Electric string states 

Reference [13] considers a class of states in Z2 lattice 
gauge theory. This theory has just two irreducible rep- 
resentations: a trivial representation and an alternating 
representation. The associated spin network states are 
"electric string states" where the two representations are 
interpreted as the presence or absence of electric strings 
along the edges. The states considered are of the form 

H = 1^E^"^"''^I^) (34) 
s 

where a is a real parameter, L{S) is the total length of 
electric string (i.e. the number of alternating representa- 
tions), and TV is a normalization factor. 

For such a state, we now show that the entangle- 
ment entropy is given entirely by the Shannon entropy 
of the representations on the boundary [the first term 
in Eq. ([32)) ]. Since the gauge group is Abelian, the sec- 
ond term in Eq. (1321) vanishes. Now consider fixing the 
set of representations on the boundary, Rg. The to- 
tal length of electric strings is the sum of the strings 



in A, the strings in B, and those crossing the bound- 
ary: L{S) = L{Sa) + L{Sb) + L{dA). For a fixed set of 
boundary representations, the reduced density matrix is 

PAiRa) cx e-^^(^-)-^^(^-M5A)(^;i] (35) 

= (36) 

where 

Sa 

This is a pure state, so S{pA{Rd)) = 0. 

Here we rederive the result that for the states in 
Eq. p4l) the entropy is just the Shannon entropy of the 
string end points. The entropy for this special class of 
states was originally derived in Ref. [13] , and on that ba- 
sis it was conjectured that the Shannon entropy of the 
string end points is a good approximation to the full en- 
tanglement entropy of the ground state in lattice gauge 
theory. The advantage of the Shannon entropy over the 
full entropy is that it depends only on the probability 
distribution of representations on the boundary, and is 
more easily computed in computer simulations. 

Here we have proven that the Shannon entropy is a 
lower bound to the entropy that appears generically for 
all states, and not just states of the special form p4l) . 
Moreover we have characterized precisely the difference 
between the Shannon entropy and the full entanglement 
entropy. The fact that Ref. fl^ finds good agreement 
between the Shannon entropy and the full entropy for 
the true ground state indicates that the third term in 
Eq. ([32)) is subleading for this state. 

For a non- Abelian gauge theory we can improve on the 
Shannon entropy as an approximation of the full entropy 
by including the log-dimension term [the second term of 
Eq. ([32)) ]. This term also depends only on the distri- 
bution of boundary representations and should therefore 
also be easy to compute in computer simulations. Refer- 
ence [13] also noted the similarity of the Shannon ent ropy 
to the log-dimension term that appeared in Ref. [21|. 
Here we have shown that these two terms are distinct 
contributions to the entanglement entropy. 

B. Topological phase ground state 

In Ref. [^, the limit a — > was considered, in which 
the state approaches a superposition of spin network 
states in which every spin network has an equal ampli- 
tude. Every configuration of string end points on the 
boundary of a region is equally probable, but gauge in- 
variance requires the total number of string end points 
crossing each connected component to be even. For a 
region with n boundary edges and whose boundary has 
k connected components, the entropy is 

S = {n-k)ln2. (38) 
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For a macroscopic region, n becomes large while k stays 
constant so the entropy is approximately extensive on the 
boundary, with small nonextensive corrections. 

The deviation from extensivity of the entanglement en- 
tropy is captured by the topological entanglement en- 
tropy. Given a pair of regions A and B,^ the topological 
entanglement entropy is (following Ref. [S], but closely 
related to the definition in Ref. [7]) 

S-t^p = S{A) + S{B) - S{A UB)- S{A n B). (39) 

Note that terms proportional to the volume or to the sur- 
face area [such as the term proportional to n in Eq. ([55]) ] 
do not contribute to the topological entanglement en- 
tropy. This is because the volume and surface area obey 
the inclusion-exclusion principle, 

f{A) + f{B)^f{ADB) + f{AnB). (40) 

where / is a function measuring either the volume or 
the surface area^. However the A;-dependent term does 
contribute to the topological entanglement entropy. If 
one considers a set of regions A, B as in Ref. such 
that A and B are each topologically disks, A fl B has 
two connected components, and AU B is topologically 
an annulus, we find a topological entanglement entropy 
of 

^top = 21n2 (41) 
in agreement with the result of Ref. Q . 

C. Strong coupling limit 

We now consider the entanglement entropy of the 
ground state of the SU(2) Kogut-Susskind Hamiltonian 
[24| in the limit of strong coupling, g 3> 1. Consider a hy- 
percubic lattice in dimension d> 2. The Kogut-Susskind 
Hamiltonian is a sum of electric and magnetic parts, 

H = He + Hb. (42) 

We will work with a rescaled version of this Hamiltonian, 
but the ground state and therefore its entanglement en- 
tropy are not sensitive to this rescaling. The electric part 
is diagonal in the spin network basis, and is given by 

He\S)=Y.3i{3i + 1)\S) (43) 
leL 

where ji is the spin of the representation r/ . The state of 
lowest energy for He is the spin network state in which 



^ We are now allowing B to be an arbitrary region, not necessarily 
the complement of A. 

^ Note that the Euler characteristic also obeys the inclusion- 
exclusion principle. This means that (counterintuitively) terms 
in the entanglement entropy proportional to the Euler character- 
istic do not contribute to the topological entanglement entropy. 



all edges are in the j = representation. We will denote 
this state by jO). 

The magnetic part of the Hamiltonian is not diago- 
nal in the spin network basis, but can be expressed as a 
functional of the holonomies, 

iJs = 3A ^ [tr(un) + h.c] (44) 

□ 

where □ is the set of all plaquettes (closed loops con- 
taining exactly four links), and ii{u\j) is the associated 
Wilson loop operator in the fundamental representation 
j = i. The parameter A is related to the gauge coupling 
g hy X ^ g^*. The operator tr(un) acts on the trivial 
spin network as 

tr(un) |0) = !□) (45) 

where ]□) is the spin network state in which each edge 
around the plaquette □ is assigned the j — ^ represen- 
tation, and all other edges are assigned the trivial repre- 
sentation. 

For strong coupling g 3> 1, so A <C 1 and we can use 
perturbation theory to calculate the ground state, treat- 
ing Hb as a perturbation of He- We will be interested 
in computing the entropy to order A^, so we compute the 
ground state to order A^: 

\n) = (l-iiVnA2)|0) + A5]|n) 

□ 

+x' I ^ !□□') + ^ cm |m) + 5] !□) ] 

\a.a' m □ / 

+0(A^). (46) 

Here iVp is the total number of plaquettes in the lattice, 
ensuring that the state is normalized to order A^. The 
state !□□') denotes the spin network state of two nonin- 
tersecting single-plaquette Wilson loops around the pla- 
quettes □ and □'. The state !□□) denotes a spin network 
state with support on two intersecting plaquettes with 
outer links in the j = \ and an intermediate link with 
j = 1, and the state ]□□) is a spin network of a single 
loop encircling two plaquettes in the j = \ representa- 
tion. The numbers cqq and are constants of order 
unity that are irrelevant for the entanglement entropy. 

To describe the way different spin networks intersect 
the region A, we will write □ G A, □□ G A to indicate 
spin networks that lie entirely in region A. We can divide 
the single-plaquette spin networks into those within A^ 
those within _B, and those intersecting the boundary. The 
numbers of plaquettes of each type are given by 7Vn(A), 
Na{B), and NQ{d), respectively, with 

Na = Na (A) + N^iB) + (d) . (47) 

To calculate N\j{d), we note that a single-plaquette loop, 
if it intersects the boundary at all must intersect an even 
number of times. We will assume that the region A is 
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chosen so that single-plaquette loops can intersect either 
twice or not at ah. Let n be the number of boundary 
Unks of region A. To count the number of ways a single 
plaquette can intersect the boundary, we fix one of the 
links intersecting the boundary, and after doing so there 
are 2{d— 1) different orientations the plaquette can take. 
This overcounts by a factor of 2, since the loop intersects 
the boundary twice, and so there are 

Na{d)=n{d-l) (48) 

ways a single plaquette can intersect the boundary. 
We now compute the entanglement entropy of the state 
by calculating each term of Eq. (|3D) in turn. To 
find the probability distribution of representations on the 
boundary we note that the probability of a two-plaquette 
state is 0{X'^) and therefore negligible. Thus the only 
states contributing to this distribution are the trivial 
spin network and the single-plaquette spin networks. The 
number of different possible sets of representations Rq is 
iVn(9) and each has probability A^, with the probability 
of having no intersections given by 1 — A^q(9)A^. The 
entropy of this probability distribution is 

H{p{Rg)) = n{d - I) {- hi + 1) + 0{X^) (49) 

and since each single-plaquette spin network intersects in 
two j = \ links, the second term of Eg. [5^ is 

(ln(2j, + I)) ^n{d^ l)X^2\n2. (50) 

leLo 

We now consider the entropy of the density matrices 
PAiRd)- Since we are taking an expectation value, we 
only need to consider sets of boundary representations 
with probability of order or larger. This means ei- 
ther there is no intersection with the boundary, or a sin- 
gle plaquette □ intersecting the boundary. In the latter 
case, the only matrix element of pA compatible with the 
assignment of representations to the boundary and prob- 
ability at least order A^ is A^ Thus pA{Rd) is a 
pure state to order A^ and so contributes no entropy. 

In the case where there is no plaquette intersecting the 
boundary, we need to know the state pA{Rd) to order A^. 
A short calculation shows that 

PA(i?9) = |^)(^|+0(A3) (51) 

where 

1^) = {l-\Nu{A)X^)\0) + xY,n (52) 

+x^{ ^ |m')+Ecm|m)+E^^I^)) • 

Since this is a pure state, its entropy is zero to order A^. 

Combining the terms in the previous paragraphs the 
entanglement entropy at first nonvanishing order in the 
strong coupling expansion is 

S = n{d- l)X^{- In A^ + 1 + 2 In 2) + 0{X^). (53) 



It is extensive in the boundary area (proportional to n). 
The entropy is also proportional to {d — 1), which is the 
number of polarizations of the gauge field. This factor is 
to be expected for weak coupling, where free field theory 
is a good approximation. It is not clear why this factor 
should appear also at strong coupling. 



IV. CONCLUSION 

We have given a formula for the entanglement entropy 
of an arbitrary state in lattice gauge theory as a sum 
of three terms [Eq. ((32)) ]. Two of these terms are local 
to the boundary and have appeared before in the liter- 
ature [13, the other captures nonlocal correlations 
between bulk degrees of freedom. Our result extends the 
result of Ref. |12| , which proposed that the Shannon en- 
tropy of the boundary representations [the first term of 
Eq. (1321) ] is an approximation to the entanglement en- 
tropy that depends only on the statistics of boundary 
observables. Our results prove that the Shannon entropy 
is a lower bound, and we give an improvement of this 
lower bound for non-Abelian gauge theories [the second 
term of Eq. ([32])] that also depends only on the statistics 
of boundary observables. Moreover, a precise expression 
is given for the difference between the local part of the 
entropy and the full entropy [the third term of Eq. (15^ ] . 

We have verified several results for entanglement en- 
tropy of specific states that appeared already in the lit- 
erature, and considered also the entanglement entropy of 
the ground state of the Kogut-Susskind Hamiltonian for 
SU(2) lattice gauge theory to first nonvanishing order 
in the strong coupling expansion. While at this lead- 
ing order only the local terms contribute to the entropy, 
at higher order all terms will contribute. This agrees 
with field theory calculations of the entropy, where the 
entropy density is found to diverge as the horizon is ap- 
proached. While we expect the dominant contribution 
to entanglement entropy to come from states localized 
near the boundary, there should be a finite contribution 
from correlations at a distance of more than one lattice 
spacing. 

It is tempting to speculate on the relation between the 
local terms in the entanglement entropy and the contact 
term found in Ref. T^l j since both appear to be unique to 
gauge theories and both are associated with observables 
localized on the boundary. However our result cannot 
explain the negative coefficient associated with the con- 
tact term, as the local part of the entanglement entropy 
is manifestly positive. There remains the intriguing pos- 
sibility that the calculation of Ref. fl^ corresponds to a 
different definition of the entropy associated with a re- 
gion of space than the one considered here, a possibility 
we leave for future work. 
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